Radiative jet energy loss in a three-dimensional hydro dynamical medium and high px 
azimuthal asymmetry of 7r suppression at mid and forward rapidity at RHIC 
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The nuclear modification factor Raa for tvo production in Au+Au collisions at yfs = 200 AGeV 
is calculated, and studied at high transverse momenta pr- The soft thermalized nuclear medium is 
described within the framework of relativistic ideal three-dimensional hydrodynamics. The energy 
loss of partonic jets is evaluated in the context of gluon bremsstrahlung in the thermalized partonic 
matter. We provide a systematic analysis of the azimuthal asymmetry of 7ro suppression at high pr 
in central and non-central collisions, at mid and forward rapidity. The determination of Raa as a 
function of pr, at different azimuthal angles, and different rapidities makes for a stringent test of our 
theoretical understanding of jet energy loss over a variety of in-medium path lengths, temperatures 
and initial partonic jet energies. This lays the groundwork for a tomography of the nuclear medium. 



I. INTRODUCTION 

Experiments at the Relativistic Heavy Ion Collider 
(RHIC) have shown that high pt hadrons in central 
A+A collisions arc significantly suppressed in comparison 
with those in binary P+p interactions, scaled to nucleus- 
nucleus collisions l2J . This phenomenon is commonly 
attributed to the fact that partonic jets produced in the 
early pre-equilibrium stage of the collisions interact with 
the hot and dense nuclear medium created in those col- 
lisions and loose energy in the process. This is referred 
to as jet-quenching Q. These lower energy partonic jets 
traverse the medium and will eventually fragment into 
hadrons which arc observed in the detectors. 

Theoretical formalisms have been elaborated to de- 
scribe the energy loss following the gluon bremsstrahlung 
experienced by the color charges in the medium: we 
mention the work by Baier-Dokshitzer-Mueller-Peigne- 
Schiff (BDMPS) Q, Gyulassy-Levai-Vitev (GLV) @, 
Kovner- Wiedemann (KW) Zakharov Majumder- 
Wang-Wang (Higher Twist) [3], and Arnold-Moore- Yaffe 
(AMY) [|. 

Recent studies [ill EH [H| indicate that additional col- 
lisional energy loss of light partons might be substan- 
tial as well. However, a consistent treatment including 
(possibly destructive) interference is to be developed [l3[ . 
In this article we will restrict ourselves to the calcula- 
tion of energy loss of the hard partons induced by gluon 
bremsstrahlung in the deconfincd phase. 

Jet quenching can be experimentally quantified by 
measurements of various quantities as e.g. the nuclear 
modification factor Raa, the elliptic flow V2 at high px, 
and high px hadron correlations. While considerable the- 
oretical effort has been deployed to develop and improve 
our understanding of modifications of jets in the nuclear 
medium, early jet quenching calculations often relied on 
an elementary description of the soft medium in their 
description of data. In most works the jets traverse a 
simple density distribution which varies with time un- 



constrained by the bulk observables. with or without a 
Bjorken expansion, see e.g. [3, [l5||. Similarly, calcu- 
lations estimating the effects of three dimensional (3D) 
expansion on Raa have treated the energy loss of jets in 
a simplified fashion [l6j . 

In [l7| a parameterized non-Bjorken fireball evolution 
that accounts for several measured observables connected 
with bulk properties of the matter created at RHIC [HI 
was applied to study the effect of flow on energy loss 
(in the BDMPS formalism, according to the prescription 
outlined in [ijl]). This study aimed for a sophisticated 
description of energy loss as well as for the medium evo- 
lution, but it was restricted to the calculation of Raa in 
central collisions. Later this observable was also studied 
in a 2D hydrodynamical evolution model [2(j. Recently, 
a 3D hydrodynamical evolution calculation 2l| of the 
expanding medium in central and non-central collisions 
was employed in detailed studies of jet energy loss as 
predicted in the BDMPS formalism 0, [23[ and in the 
higher twist formalism [24|. 

The present work contributes to this effort of under- 
standing the physics of jet quenching by applying the 
Arnold, Moore, and Yaffe (AMY) formalism Q for gluon 
bremsstrahlung to calculate the jet energy loss in the 
thermal partonic medium in central and non-central col- 
lision as inferred from 3D relativistic hydrodynamics [2l[ . 
We present a calculation of Raa as a function of trans- 
verse momentum (and the azimuth) in central and non- 
central collisions and also study the rapidity dependence 
of this quantity. 

While Raa as measured in central collisions alone is 
not suited to distinguish in detail between different theo- 
retical conjectures about jet energy loss (2f|, the combi- 
nation with additional measurements of Raa versus re- 
action plane in non-central collisions [26| and at finite 
rapidity provides further valuable tomographic informa- 
tion. Additional tomographic observables are high p t 
triggered correlation measurements, see e.g. [20l. [271 [28j . 

The paper is organized as follows, we first briefly re- 
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view the 3D hydrodynamical description of the medium 
in Sec. II. We then discuss in Sec. Ill how the initial mo- 
mentum distributions of jets and their time-evolution in 
the thermal medium (which incorporates the energy loss 
process in the AMY formalism) as well as the fragmenta- 
tion of the final jets into pions are calculated. Numerical 
results are presented for Raa at mid and forward rapid- 
ity in Sec. IV together with a comparison to data where 
already available. Finally, Sec. V contains our conclu- 
sions. 



number density), respectively. W(x,y;b) is given by 
a combination of wounded nucleon model and binary 
collision model [32| and H(rf) is given by H(rf) = 
eiq ? [-(\r 1 \- Vo ) 2 /(2^)-e(\r ] \- Vo )]. 

The initial conditions have been chosen such that a 
successful description of the soft sector at RHIC (elliptic 
flow, pseudo-rapidity distributions and low-px single par- 
ticle spectra) is achieved. For further details, especially 
also a discussion of the EoS which is employed, we refer 
the reader to [2l| . 



II. 3D HYDRODYNAMICAL MEDIUM 



III. JET EVOLUTION AND FRAGMENTATION 



The behavior related to the bulk properties of the high- 
density phase in heavy-ion collisions at RHIC is well de- 
scribed by Relativistic Fluid Dynamics (RFD, see e.g. 
[29l . [30L HH), while this description is not applicable in 
the late dilute stages of the collisions in which the mean 
free path of hadrons is large on the typical scales of the 
system. 

In the present paper we use a fully 3D hydrodynami- 
cal model for the description of RHIC physics [2l[ which 
solves the relativistic hydrodynamical equation 



d,T^ = 0, 



(1) 



where T^ v is the energy momentum tensor which can be 
expressed as 



T" v = (e + p)U»U v -pg"". 



(2) 



Here e, p, U and <? Ml/ are energy density, pressure, four 
velocity and metric tensor, respectively. Furthermore 
baryon number n B conservation is imposed as a con- 
straint 



d^n B (T,ii)Un = 0, 



(3) 



and the resulting set of partial differential equations is 
closed by specifying an equation of state (EoS): e = e(p). 
Our particular RFD calculation utilizes a Lagrangian 
mesh and the coordinates (r, x, y, rj) with the longitu- 
dinal proper time r = \/t 2 — z 2 and space-time rapidity 
rj = |ln[(£ + z)/(t — z)] in order to optimize the calcu- 
lation for the ultra-rclativistic regime of heavy collisions 
at RHIC. Once an initial condition has been specified 
RFD in the ideal fluid approximation (i. e. neglecting 
off-equilibrium effects) allows a calculation of single soft 
matter properties at RHIC, especially collective flow and 
particle spectra. 

We assume early thcrmalization with subsequent hy- 
drodynamical expansion at tq = 0.6 fm/c. The ini- 
tial conditions, namely initial energy density and baryon 
number density arc parameterized by 



n B (x,y,rj) 



W(x,y; b)H(r)), 
n Bmax W(x,y;b)H(r/), 



(4) 



where b and e max (n Bmax ) arc the impact parame- 
ter and the maximum value of energy density (baryon 



In this section we present the techniques used to cal- 
culate the initial jet production in the early stage of the 
collisions, the subsequent propagation through the hot 
and dense medium, and final hadronization in the vac- 
uum. We exclusively focus on the hadrons in the high 
Pt region in which fragmentation is the dominant mech- 
anism for the production of hadrons. For softer hadrons 
(below pt ~ 7 GcV/c) other mechanisms, such as the re- 
combination of partons become of increasing significance 

M- 

The initial jet density distribution VabO 3 , r_L) at the 
transverse position rj_ and in A+B collisions with impact 
parameter b is given by 



VA B (b,r ± ) 



T A (f±+b/2)T B (f ± -b/2) 
Tab {b) 



(5) 



Here we use a Woods-Saxon form for the nuclear den- 
sity function, p(r±,z) = po/[l + cxp( 1 ^)], to evaluate 
the nuclear thickness function Ta{t±) = J dzpA{r±, z) 
and the overlap function of two nuclei TA B (b) = 
/ d 2 r±TA(r±)T B (f± + b). The values of the parameters 
R = 6.38 fm and d = 0.535 fm are taken from (34|. 

The initial momentum distribution dN AB (b) / d 2 p T dy\i 
of jets is computed from pQCD in the factorization for- 
malism, 



dN AB (b) 



d 2 p T dy 



TAB(b) ^ J dx a G a /A(x a ,Q)G b/B (x h ,Q) 
r dcr a +b^j+d 



1 



abd 



f 2x a - xte y 



-K- 



dt 



(6) 



In the above equation, G a /A{x a ,Q) is the distribution 
function of parton a with momentum fraction x a in the 
nucleus A at factorization scale Q, taken from CTEQ5 
35 1 including nuclear shadowing effects from EKS98 
36j . The index j represents one of the partonic species 

(j = <lT<j>g): an d x T = 2p T j \Jsnn, where ^/smn is the 
center of mass energy. The distribution da/dt is the 
leading order QCD differential cross section, and the K- 
factor accounts for NLO effects and is taken to be con- 
stant in our calculation as it is almost Pt independent 
[iiTl . liiR .'!!) . The initial Cronin effect is neglected in our 
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calculation since the nuclear modification factor of neu- 
tral pions from d+Au collisions measured by PHENIX is 
consistent with 1 within systematic errors [40]. 

The evolution of a jet momentum distribution 
Pj (p, t) = dNj (p, t) I dpdy (essentially the probability of 
finding a jet with energy p at time t) in the medium is 
obtained in the AMY formalism by sol ving a set of cou- 
pled rate equations (for details see [4j], IHf) , which have 
the following generic form, 



dt 



= £/* 

ab 



P a (p + k,t) 
drj b ( P ,k) 

dkdt 



dT%(p + k,k) 



dkdt 



(7) 



where dT 3 ab (p, k) / dkdt is the transition rate for the par- 
tonic process j — ► a+b. Wc point out that the calculation 
includes not only the emission but also the absorption of 
thermal partons as the k integral in Eq. ( |3 ra nges from 
— oo to oo. The transition rate is given by [4ll. |42| 



dT(p,k) 
dkdt 



Cs9l 



1 



1 



Wirp" 7 1 ± e- k / T 1 ± e-(P- k )/ T 



x < 



x 3 (l — x) 2 



N { 



2 +(l-*) 2 

x 2 (l-x) 2 
l+x i + (l~x) i 
x 3 (1 — x) 3 




d 2 h 

(2tt)' 



2h-ReF(h,p, k) 



(8) 



Here C s is the quadratic Casimir relevant for the pro- 
cess, and x = k/p is the momentum fraction of the gluon 
(or the quark, for the case g — > qq). h = p X k deter- 
mines how non-collinear the final state is; it is treated 
as parametrically 0(g s T 2 ) and therefore small compared 
to p ■ k. Therefore it can be taken as a two-dimensional 
vector in transverse space. F(h, P, k) is the solution of 
the following integral equation [4l|, [42[ : 



2h 



i5E(h,p, k)F{h)+g 2 s 



d 2 q± 



C(qi 



(2n) 2 

:{(C S - C A /2)[F(h) - F(h-kq±)) 

-(C A /2)[F(h)-F(h+pq ± )) 
-(C A /2)[F(h) - F(h-(p-k)q ± )}} 



(9) 



Here 5E(h,p, k) is the energy difference between the final 
and the initial states: 



6E(h,p, k) 



h 2 



2pk(p—k) 



2k 



" L P -k 
2(p-k) 



2p 



and to 2 are the medium induced thermal masses. Also, 
C(q±) is the differential rate to exchange transverse (to 
the parton) momentum q±. In a hot thermal medium, 



its value at leading order in a s is [4c 



6 



(2N C +N { ). (11) 



For the case of g — > qq, (C s — C\/2) should appear as 
the prefactor on the term containing F(h — pq±_) rather 
than F(h — kq±). 

The strength of the transition rate in pQCD is con- 
trolled by the strong coupling constant a s , temperature 
T and the flow parameter /3 (the velocity of the thermal 
medium) relative to the jet's path. In a 3D expanding 
medium, the transition rate is first evaluated in the lo- 
cal frame of the thermal medium, then boosted into the 
laboratory frame, 



dT(p, k) 



dkdt 



lab 



dkndt 



O aL 



local 



where k Q = k(l - Vj ■ f3)/^Jl-f3 2 and t = t^X-0 2 
are momentum and the proper time in the local frame, 
and Vj is the velocity of the jet. As jets propagate in the 
medium, the temperature and the flow parameter depend 
on the time and the positions of jets, and the 3D hydrody- 
namical calculation pH is utilized to determine the tem- 
perature and flow profiles. The energy- loss mechanism 
is applied at tq = 0.6 fm/c, when the medium reaches 
thermal equilibrium, and turned off when the medium 
reaches the hadronic phase. 

The final hadron spectrum dN^ B {b) / d 2 prdy at high 
Pt is obtained by the fragmentation of jets in the vacuum 
after their passing through the 3D expanding medium, 



d 2 prdy 



J d 2 r 1 _V AB {b,r 1 _) J -^-D h/j ( Zj ,Q F ) 

i 3 



dN\ B (b,r ± ) 



d 2 pLdy 



(13) 



where dN 3 AB (b, f±_)/d 2 p 3 T dy\f is the final momentum dis- 
tribution of the jet initially created at transverse posi- 
tion rj_ after passing through the medium. The frag- 
mentation function D h /j(zj,Qp) gives the average mul- 
tiplicity of the hadron h with the momentum fraction 
zj = pt/p j t produced from a jet j at a scale Qf, taken 
from KKP parametrization [44( . The fractorization scale 
Q = pip and fragmentation scale Qf = Pt are set as in 
[371 ] where the X-factor is found to be 2.8. We use these 
values throughout the present study. This nicely repro- 
duces the experimentally measured 7To yield at mid and 
forward rapidity in p+p collisions at yfs = 200 GcV, as 
shown in Fig. [T] and Fig. [2] It can be clearly seen that 
replacing the CTEQ5 parton distribution functions by 
MRST01 [i^ yields essentially the same result for the 
inclusive ttq production in p+p collisions. We point out 
that the presence of a nuclear medium might in principle 
alter these scales but we postpone a detailed study of this 
possibility to future research. 



4 



10" 



10' 



L0 



S to" 4 



<•£ io- 6 



10" 



10 



• PHENIX, Preliminary 
CTEQ5, KKP, Q = p T J , Q F = p T , K = 



MRST01.KKP, Q = 



, K = 2.8 



10 

p T (GeV/c) 



15 



20 



FIG. 1: (Color online) The inclusive cross section for ttq, pro- 
duction versus no transverse momentum at mid-rapidity in 
pp collisions at y/s = 200 GeV, compared with PHENIX data 



200 AGcV. We only show results for the nuclear mod- 
ification factor Raa for neutral pions, as results for 
charged hadrons (including contributions from charged 
pions, kaons and protons) are qualitatively similar. 
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FIG. 2: (Color online) The inclusive cross section for ttq pro- 
duction versus no energy at forward rapidity in pp collisions 
at yfs = 200 GeV. Data points are taken from STAR [13]. 



The nuclear modification factor Raa is defined as the 
ratio of the hadron yield in A+A collisions to that in p+p 
interactions scaled by the number of binary collisions 



RAA(b,PT,y) 



1 



dNX A (b)/d 2 p T dy 



Nco U (b) dN^/d?p T dy 



(14) 



IV. RESULTS 

In Fig. [3J we present the calculation of the nuclear 
modification factor Raa for neutral pions measured at 
midrapidity for two different impact parameters b = 
2.4 fm and b = 7.5 fm, compared with (preliminary) 
PHENIX data for most central (0 — 5%) and midcen- 
tral (20 — 30%) collisions (26|. All results presented 
throughout the paper are for Au+Au collisions at = 



FIG. 3: (Color online) The neutral pion Raa at midrapidity 
in most central (upper panel) and midperipheral (lower panel) 
Au+Au collisions compared with PHENIX data. 

Once the temperature evolution is fixed by the initial 
conditions and subsequent 3D hydrodynamical expan- 
sion, the strong coupling constant a s is the only quantity 
which is not uniquely determined in the model. The value 
of a s is a direct measure of the interaction strength and 
is adjusted in such a way that the experimental data 
in the most central collisions is described. The same 
value, a s = 0.33, is used in peripheral collisions. Treat- 
ing a s as T-independent from early thermalization on 
down to the phase transition temperature is a simplifica- 
tion and corresponds to the assumption that the decon- 
fined phase of the medium formed in Au+Au collisions at 
200 AGeV at RHIC can be characterized by one average 
effective coupling. We point out that a treatment in the 
AMY formalism only considers energy loss in the par- 
tonic phase, hadronic energy loss is not included in the 
present study [49| . We have verified that choosing dif- 
ferent constant values of a s does not influence the shape 
of Raa as a function of pr significantly while only the 
overall normalization is affected. We point out that al- 
though a 8 < 1, g s = \/4ira s is actually larger than 1. In 
that sense our study does not contradict the finding in 
[22| that a stronger quenching power of the medium has 
to be assumed than if a fully perturbative treatment of 
let quenching in the quark gluon plasma is employed. (In 
22l ] uncertainties in the selection of the strong coupling 
and possible non-perturbative effects were parameterized 
by a factor K in q = 2Ke 3 ^ 4 . K w 3.6 was adjusted to 
give a good description of Raa in central collisions at 
mid-rapidity.) 

In Fig. [31 Raa a t midrapidity is averaged over the az- 
imuth </>. More tomographic capabilities can be achieved 
if one studies Raa at midrapidity in non- central colli- 
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sions not only as a function of px averaged over <fi but 
also as a function of the azimuth cf> [2(| . 

The reason is that the initial geometric asymmetry in 
non-central collisions leaves its imprint on the 3D hydro- 
dynamical evolution and initial jets experience different 
energy loss (depending on where they are produced in the 
medium and in which direction they are emitted) owing 
to the different local properties of the nuclear medium 
with which they interact. In the AMY formalism the im- 
portant input from the evolution is the temperature in 
the rest frame of the local fluid that the jet experiences 
(and to a lesser extent the flow profile of the medium, as 
discussed later). 



t = 0.6 fm/c, inner 

t = 0.6 fm/c, outer 

— t = 3.6 fm/c, inner 

— t = 3.6 fm/c, outer 




4 5 6 
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FIG. 4: (Color online) The inner and outer boundaries for 
T = T c in the transverse plane at two different proper times, 
b = 7.5 fm. 
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times of the evolution. They represent the inner and 
outer boundaries of the mixed phase during the evolu- 
tion. The geometric asymmetry of the temperature pro- 
file can be clearly seen from the plot. Both boundaries 
move towards the center and the inner boundary moves 
faster than the outer boundary. It is useful to define the 
emission in plane (cf> = 0) versus out of plane (0 = 7r/2). 
We point out that the ratios of the boundary positions in 
plane to those out of plane are almost constant in proper 
time and have almost the same values for the inner and 
outer boundaries: ~ 0.8. 

Fig.[5]shows the temperature observed by a jet travers- 
ing this medium. The jet is assumed to be created at 
position (ro,4>o) by a hard scattering at early times in 
the heavy-ion collision. As it propagates through the 
medium, the surrounding environment will change from 
the QGP phase to the mixed phase, then to the hadronic 
phase and will eventually freeze-out. We plot the tem- 
perature evolution experienced by jets that are created in 
a symmetric position (<fio = 7r/4) relative to in-plane and 
out-of-plane and illustrate the geometrical asymmetry of 
the medium. We compare jets starting at the origin and 
those at r = 3 fm. 

Jets that propagate out of plane will pass the mixed 
phase and the hadronic phase at later proper time than 
those traversing in plane and will interact with the de- 
confined and mixed phase of the medium longer. If the 
jets have identical initial energy, the energy loss experi- 
enced by the jets propagating out of plane will therefore 
be larger than in plane. 
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FIG. 6: (Color online) The neutral pion Raa at midrapidity 
for emissions in plane and out of plane as a function of pr for 
different impact parameters. 



FIG. 5: (Color online) The time evolution of the temperature 
seen by a jet initially created at (ro, 4>o) moving in plane and 
out of plane through the medium, b — 7.5 fm. 

To illustrate the geometrical asymmetry we show in 
Fig. [3] isotherms for T = T c in the transverse plane for an 
impact parameter of b = 7.5 fm at two different proper 



This behavior is reflected in Raa as a function of pt for 
emissions in plane and out of plane in Fig. [S] for different 
impact parameters. While there is very small difference 
for Raa between the two planes in central collisions, a 
much larger difference for midccntral collisions (about 
13% for b = 7.5 fm) is predicted, as can been seen from 
the ratio of Raa for emission out of plane to that in plane 
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FIG. 7: (Color online) The ratio of the neutral pion Raa 
at midrapidity for emissions in plane and out of plane as a 
function of pr for different impact parameters. 



FIG. 9: (Color online) Comparing the neutral pion Raa at 
midrapidity with and without flow for emissions in plane and 
out of plane as a function of pr, b — 7.5 fm. 



as shown in Fig. [7J 
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FIG. 8: (Color online) The neutral pion Raa at midrapidity 
as a function of the azimuthal angle <f> of the pion for different 
pr, b = 7.5 fm. 

As a further tomographic quantity, one can also 
study Raa for non-central collisions as a function of 
the azimuthal angle <f> for different pt, see Fig. [5] A 
monotonous decrease of Raa for emissions from in plane 
to out of plane, reflects (an average of) the asymmetric 
temperature (and flow) profiles experienced by the jets 
while they traverse the medium. 

In a 3D expanding medium, there is also considerable 
collective flow being built up during the evolution. This 
can affect the energy loss of jets and may to some de- 
gree influence the asymmetry in the final pion spectrum. 
To quantify this effect, we use the same 3D hydro tem- 
perature profile, but disregard the transverse flow. We 
compare the case with flow to one where the velocity ef- 
fect is disregarded, namely /3 = is enforced by hand in 
Eq. (fT2")) (only for illustration purposes). This treatment 



can give an estimate on how collective flow (not the tem- 
perature of the medium) influences the jet energy loss in 
the evolution. As is shown in Fig. [9l How effects only 
slightly increase the quenching power of the medium in 
the AMY-formalism. It is emphasized that for a realistic 
hydrodynamical calculation, the overall temperature of 
the medium would drop not as fast if collective flow was 
switched off and the medium itself would expand more 
slowly in this case. 

We point out that a further interesting quantity is Raa 
for neutral pions as a function of px at different cen- 
tralitics and away from midrapidity. The formalism as 
outlined in Section II. can be straightforwardly extended 
to treat this case. Caveats are that only moderate de- 
viations from midrapidity can be allowed, because the 
nuclear parton distribution functions can be less exactly 
determined in the relevant region (36j and the assumption 
of a thermalized medium essential for a hydrodynamical 
treatment is no longer fulfilled far away from midrapid- 
ity. We therefore restrict our study to rapidities close to 
midrapidity (maximum forward rapidity y = 2). 

At finite rapidity y the energy of a highly-relativistic 
jet with a transverse momentum px is given by E = 
Pt cosh y. The pions at a fixed px have more energy and 
are the fragments of higher energetic partons than the 
corresponding midrapidity pions. The initial jet distri- 
bution of quarks and anti-quarks is shown in Fig. [TO] for 
different rapidities, compare Eq. ©. Note that the kine- 
matical cut off at E = t/sJjn '/2 = 100 GeV is reached at 
lower pt for finite y. 

In Fig. [TT]we show Raa as a function of pt for central 
collisions (0 — 5%, b = 2.4 fm) at mid and forward ra- 
pidity. It is interesting to notice that Raa behaves quite 
differently as a function of px at y = 2 than at y = 0. 
This is not only due to the different temperature profiles 
of the hydrodynamical medium at forward rapidity but 
also strongly influenced by the different initial jet distri- 
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FIG. 10: (Color online) The jet (quark + anti-quark) trans- 
verse momentum distribution at different rapidities, b — 
2.4 fm. 
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FIG. 12: (Color online) Comparing neutral pion Raa with 
and without nuclear shadowing effect at different rapidities, 
b = 7.5 fm. 



jet distribution according to Eq. ([6]) at high transverse 
momentum which decreases faster than an overall power 
law, see Fig. [TO] 
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To provide additional insight, we studied the same 
quantity averaged over <f> for midcentral collisions with 
an impact parameter of b = 7.5 fm with and without 
nuclear shadowing effects taken into account in the par- 
ton distribution functions utilized in Eq. @. Results 
are shown in Fig. [T2J It is interesting to notice that 
Raa is not monotonously increasing as a function of px- 
The midrapidity Raa is decreasing above ~ 18 GeV/c 
(with nuclear shadowing), the turning point for y = 1 is 
at ~ 9 GeV/c (with nuclear shadowing). The values of 
Raa at y = 2 decreases monotonically above ~ 6 GeV/c 
in the case without nuclear shadowing and exhibits two 
turning points if shadowing is taken into account. We 
have also found that assuming a simple power law ap- 
proximation for dN I 'd prdy distributions for all values 
of pt would lead to increased Raa at higher px (com- 
parison not shown). This demonstrates that the overall 
decrease of Raa at higher px is mainly due to the initial 



FIG. 13: (Color online) The ratio of the neutral pion Raa im- 
posing a boost-invariant expansion to Raa as calculated from 
the 3D hydrodynamical (non-Bjorken) medium b — 7.5 fm. 

One can also address the question how different the 
part of the medium is traversed by a jet which frag- 
ments into a pion at forward rapidity in comparison to 
one which fragments at midrapidity. We compare the full 
3D hydrodynamical calculation to an effective 2D boost- 
invariant approach in which the 2D hydrodynamical solu- 
tion at midrapidity is assumed to also describe the trans- 
verse profile at forward rapidity. This corresponds effec- 
tively to imposing a posteriori Bjorken expansion onto 
the non-Bjorken hydrodynamical evolution. We study 
the ratio of Raa by imposing a boost invariant expan- 
sion, and comparing with the fully 3D non-Bjorken evo- 
lution. Fig. [T3] shows a calculation at forward rapidities 
for non-central collisions with a finite impact parameter 



8 



of b = 7.5 fm. This ratio is obviously not measurable, 
but is interesting from a theoretical point of view. Its 
relatively strong deviations from 1 at y = 2 stem mainly 
from the different transverse temperature profiles at for- 
ward rapidity in the non-Bjorken evolution whereas these 
differences at y = 1 are not significant. The fact that the 
ratio is rather flat in pt (it varies only in the range of 
0.7 ± 0.05 for y = 2) indicates that the reduction of the 
quenching power of the medium in the non-Bjorken case 
compared to the boost-invariant one is similar for par- 
tons over the full range of initial jet energies probed in 
the collision. Therefore a measurement of the absolute 
normalization of Raa at midrapidity and forward rapidi- 
ties might be useful in quantifying the deviations arising 
from the simplifications made in boost invariant expan- 
sion models. 



Dependence on Nuclear Parton Distribution 
Functions 
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FIG. 14: (Color online) The neutral pion Raa at midrapidity 
in most central (upper panel) and midperipheral (lower panel) 
Au+Au collisions compared with PHENIX data. Different 
prescriptions of nuclear parton distribution functions are used 
for comparison. 

As has been pointed out earlier, see e.g. [48[, the de- 
termination of nuclear parton distribution functions (nu- 
clear PDFs) from experimental data is ambiguous. These 
uncertainties can also influence the calculation of the nu- 
clear modification factor at mid and forward rapidity at 
RHIC. We compare in this subsections results obtained 
with the nuclear parton distribution functions as deter- 
mined by NPDF04 [48[ with those that were employed 
so far in this work, namely EKS98. We checked that the 
nucleon parton distributions which NPDF04 and EKS98 
rely on, namely MRST01 and CTEQ5, respectively, lead 
to almost the same prediction of the inclusive cross sec- 
tion for 7r production in p + p collisions, cmp. Figs. [T] 
andO This should be expected since the determination 
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FIG. 15: (Color online) Comparing neutral pion Raa at dif- 
ferent rapidities using different descriptions of nuclear parton 
distribution functions, b — 2.4 fm. 



of nucleon parton distribution functions has smaller un- 
certainties than those extended to nuclei. 

Fig. [14] shows the neutral pion Raa at midrapidity in 
central and midperipheral Au+Au collisions as obtained 
with the two different nuclear PDFs. Differences due to 
the different nuclear PDFs appear especially at larger 
transverse momenta of the produced pions. The same 
holds true for Raa at forward rapidity, see Fig. [T5] for a 
comparison in central collisions. 
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FIG. 16: (Color online) The ratio (NPDF04/EKS98) for the 
initial quark + anti-quark jet distributions (upper) and nu- 
clear modification factors Raa (lower) using different descrip- 
tions of nuclear parton distribution functions, b — 2.4 fm. 

It is possible to trace these differences in Raa back to 
differences in the initial jet distributions resulting mainly 
from the different shadowing descriptions. We show in 
Fig. [16] upper panel the ratio of the initial quark and 
anti-quark jet distributions as inferred from NPDF04 to 
EKS98. This translates - after jet-energy loss and frag- 
mentation have been taken into account - into a similar 
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behavior of the ratios of the nuclear suppression factor 
Raa hi the two cases. Differences in the initial distribu- 
tion (mainly resulting from different nuclear shadowing) 
will therefore be reflected in Raa at mid and forward 
rapidity and at different centralities. Reduced sensitivity 
to differences in shadowing effects is expected if ratios 
of Raa are considered: the ratio for y = 1 (Fig. fT5|) 
is only sensitive on the 1% level to employing EKS98 
or NPDF04 (comparison not shown), the sensitivity for 
y = 2 is at most 4%. What we find in this subsection 
clearly demonstrate that that Raa is not only sensitive 
to the employed jet quenching formalism but also to nu- 
clear shadowing effects. The reason is that - even after 
energy loss and fragmentation - Raa is sensitive to the 
initial jet distribution which in turn vary within the un- 
certainties of the determination of nuclear shadowing. A 
further reduction of uncertainties in the determination of 
nuclear shadowing effects will make a more stringent test 
of jet quenching formalisms by Raa measurements. 

V. CONCLUSIONS 

In this paper, the jet energy loss was studied in the 
AMY-formalism using a 3D hydrodynamical evolution 
model that has been shown to describe the bulk proper- 
ties of matter created in heavy-ion collisions at RHIC. 

We have evaluated the nuclear suppression factor Raa 
for neutral pions in central collisions as a function of 
Pt at mid and forward rapidity and have discussed how 
the azimuthal asymmetry of the medium in non-central 
collisions allows to put stronger constraints on our un- 
derstanding of jet energy loss by gluon radiation. Since 
the jets probe different flow and temperature profiles in 
the asymmetric expansion depending on their initial po- 
sitions and emission angles, Raa is not only a function 
of Pt but also of the azimuth in those collisions. 

The measured Raa as a function of pr in central and 



(averaged over the azimuth) in non-central collisions at 
midrapidity is in good agreement with the model calcu- 
lations. We also have provided calculations of Raa as a 
function of pt and the azimuth without averaging that 
can be the basis of more stringent experimental tests once 
further data become available. We furthermore studied 
Raa as a function of pt in central and again averaged 
over the azimuth in non-central collisions at mid and for- 
ward rapidity and provided arguments that a measure- 
ment of these dependences might not only be able to 
reveal more information about the nuclear medium (as 
deviations from the assumption of boost invariance) but 
also provide a possibility to observe nuclear shadowing ef- 
fects in the initial parton distribution function indirectly 
(assuming appropriate experimental resolution). 

We emphasize that the description of Raa (as a func- 
tion of pt , the azimuthal angle and rapidity) alone is not 
enough to prove the consistency of a specific energy loss 
mechanism with data, if assumptions about the medium 
evolution can be freely adjusted. On the contrary, Raa 
will only provide stronger constraints on our theoretical 
physical conjectures about jet energy loss in the nuclear 
medium if studied in a dynamical evolution model which 
has been tested using soft observablcs. 
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